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ABSTRACT 



^ ■ A large portion of the Galaxy {£ = 123° to 164°, b = —6° to —35°), which samples 

£J2 ■ regions of the Local (Orion) spiral arm and the more distant Perseus arm, has been 

■ mapped with the Wisconsin H-Alpha Mapper (WHAM) in the [S II] A6716 and [N II] 

I A6583 lines. By comparing these data with the maps from the WHAM Ha Sky 

Q>^ I Survey, we begin an investigation of the global physical properties of the Warm Ionized 

' Medium (WIM) in the Galaxy. Several trends noticed in emission-line investigations 

• of diffuse gas in other galaxies are confirmed in the Milky Way and extended to 
much fainter emission. We find that the [S II] /Ha and [N II] /Ha ratios increase as 



absolute Ha intensities decrease. For the more distant Perseus arm emission, the 



c/2 , increase in these ratios is a strong function of Galactic latitude, b, and thus, of height, 

^ ' z, above the Galactic plane, while the [S II]/[N II] ratio is relatively independent of 

Ha intensity. Scatter in this ratio appears to be physically significant, and maps of 
^ ' [S II] /[N II] suggest regions with similar ratios are spatially correlated. The Perseus 

■ arm [S II]/[N II] ratio is systematically lower than Local emission by 10%-20%. 

With [S II] /[N II] fairly constant over a large range of Ha intensities, the increase of 
[S II] /Ha and [N II] /Ha with \z\ seems to reflect an increase in temperature. Such 
an interpretation allows us to estimate the temperature and ionization conditions in 
our large sample of observations. We find that WIM temperatures range from 6,000 
K to 9,000 K with temperature increasing from bright to faint Ha emission (low to 
high [S II] /Ha and [N II] /Ha) respectively. Changes in [S II] /[N II] appear to reflect 
changes in the local ionization conditions {e.g. the S+/S++ ratio). We also measure 
the electron scale height in the Perseus arm to be 1.0 it 0.1 kpc, confirming earlier, less 
accurate determinations. 



Subject headings: ISM: structure — ISM: atoms — HII regions — Galaxy: halo — 
ISM: general 
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1. Introduction 

Although the Wisconsin H-Alpha Mapper (WHAM) Sky Survey will present the first view of 
the distribution and kinematics of the Warm Ionized Medium (WIM), we still know little about 
the physical conditions in this gas or how they are produced. Several studies have measured 
constraints on important parameters toward specific regions of the WIM, but none has yet 
provided a large sample of measurements to help test how environmental conditions {e.g. distance 
from H II regions, distance from the Galactic plane, etc.) might affect these physical conditions. 

Reynolds (1985a) derived the first estimate of the temperature of the WIM by combining 
measured line widths of [S II] and Ha emission toward the same regions. He was able to place 
limits between 5,000 K and 20,000 K, with a mean of about 8,000 K, on the temperature of a 
sample of 21 WIM observations. Ionization conditions have been probed even less, with most of 
our current information coming from observations in the plane of the Galaxy. Reynolds (1985b) 
showed that [O III] A5007 emission is very faint in the WIM, placing limits on the amount of 
0++ (needing hiy > 35.1 eV) in the WIM. More conclusively, Reynolds & Tufte (1995) and Tufte 
(1997) have shown via observations of the He I A5876 recombination line that the fraction of He"*" 
(needing hiy > 24.6 eV) in this gas is also quite low in comparison with traditional O star H II 
regions. The recent detection of [O I] A6300 from three WIM directions (Reynolds, et al. 1998) 
places further limits on the ionizing spectrum. The low [O I] /Ha ratios seen in all three directions 
suggests that the H is mostly H"*" in the Ha emitting gas. 

The advent of WHAM provides a unique opportunity to trace some of these diagnostics over 
a larger region of the sky. This paper presents the first velocity-resolved study of [S II] and [N II] 
from the diffuse gas of our Galaxy. A region of the sky that includes the Perseus arm was chosen 
for examination since the Ha emission from the WIM in this arm is well separated in velocity from 
local WIM emission and can be traced to high Galactic latitude, allowing an exploration of the 
emission with distance, z, from the Galactic plane. Such observations can be directly compared to 
studies of edge-on galaxies, most notably those by Rand (1997, 1998) of NGC 891. The primary 
trends found in his work are borne out here and extended to even fainter emission measures. 
Specifically, [S II] /Ha and [N II] /Ha ratios increase with \z\, while [S II] /[N II] ratios remain nearly 
constant. 

A brief description of the instrument and observations are presented in The results 
from the data, including a measurement of the electron scale height in the Perseus arm and an 
investigation of the trends in line ratios, are presented in §^ In we present our interpretation 
that these line ratios provide direct measurements of the temperature and ionization state of the 
WIM. Finally, our conclusions are summarized in 



2. Observations 
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2.1. The WHAM Instrument 

Due to the extended nature of the WIM, large-aperture Fabry-Perot detection techniques have 
proven to be quite successful in studying the WIM (Reynolds et al. 1990). For WHAM, a 6-inch, 
dual-etalon Fabry-Perot spectrometer delivers the high-efficiency, high-resolution observations 
(Tufte 1997; Reynolds et al. 1998). The dual-ctalon design provides superior order rejection 
over singlc-ctalon systems to achieve a much larger free spectral range while maintaining high 
spectral resolution (12 km s~^). An all-sky siderostat housing a 0.6-meter lens guides light to 
the etalon chambers in an environmentally controlled trailer. The center wavelength of the 200 
km s~^ spectral window is selected by changing the pressure of a high index of refraction gas 
(SFg) in sealed etalon chambers. A narrow-band (~ 20 A FWHM) interference filter provides 
additional spectral isolation near the desired wavelength. The resulting spectrum, in the form of a 
Fabry-Perot ring pattern, is imaged onto a cryogenically cooled Tek 1024x1024 CCD. 

The spectrum covers a 200 km s"^ interval at 12 km s^^ spectral resolution. Coatings on 
the etalons and optics allow the placement of the 200 km s~^ wide spectral window anywhere 
between 4800 and 7300 A. In this "spectral mode," the image contains only the average spectrum 
within the one-degree diameter circular beam, avoiding confusion between spectral features and 
angular structures of the source (or stars) within the beam. WHAM also has a set of optics that 
can be placed in the post-etalon beam to provide an extremely narrow-band (adjustable between 
12-200 km s~^), 1' angular resolution image of the sky within a onc-dcgrcc field of view. Only the 
spectral mode (with one-degree angular resolution) is used for the observations presented here. 
Further details on the WHAM spectrometer and its calibration can be found in Tufte (1997). 

WHAM is currently located at Kitt Peak National Observatory. Using the WIYN messaging 
system designed by Percival (1994, 1995) and custom hardware designed by the Space Astronomy 
Lab at Wisconsin, WHAM is a completely remotely operated telescope system. All functions of 
the siderostat and spectrometer are routinely operated from our Wisconsin offices. Details on 
the hardware and software used to operate WHAM and to automate the survey can be found in 
Haffner (1999). 

2.2. The Data 

The Ha observations of the 1100 deg^ region i = 123° to 164°, b = -6° to -35° that is the 
focus of this study were obtained between 1997 July and 1997 October as a part of the WHAM 
Ha Sky Survey (Haffner 1999). Additional observations at the same survey grid points but with 
the 200 km s~^ spectral window tuned near the [S II] A6716 and [N II] A6583 lines were obtained 
in 1997 October, 1997 November, and 1998 June. The region presented here is sampled by 1360 
Ha, 1360 [S II], and 1360 [N II] spectra, each representing the integrated emission within WHAM's 
one-degree beam. The integration time for each spectrum was 30 s for Ha and 60 s for the other 
two emission lines. 
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The [S II] and [N II] observations were converted into spectra using the same steps as those 
for Ha presented in Haffner (1999). Separate white-light flat fields were created for each of 
these additional filters and used in place of the Ha flat field. There are no strong terrestrial 
components in [S II] or [N II] comparable to the geocoronal emission in the Ha spectra. For 
the [S II] observations, the weak terrestrial (Tucson lights) Ne I line at 6717.04 A was used for 
velocity calibration. For the [N II] observations, spectra of Ha and [N II] from bright emission 
nebulae were compared and the velocity frame of [N II] was adjusted so that the mean velocity 
of the nebular [N II] emission matched that of the Ha. By using the terrestrial lines in the Ha 
spectra (geocoronal Ha) and [S II] spectra (Ne I), the velocity scale in each spectrum is calibrated 
to better than 1 km s~^. Intensities were calibrated by applying a correction factor to the [S II] 
and [N II] observations (0.94 and 1.15, respectively) based on the transmission differences between 
these filters and the Ha filter. 

As in the Ha observations, the [S II] and [N II] spectra are contaminated by weak atmospheric 
lines (Haffner, Reynolds, & Tufte 1998; Haffner 1999). The brightest of these in our [S II] spectra 
is the Ne I line located +26.8 km s~^ from the geocentric zero of [S II] A6716. As noted above, it 
proves to be a useful wavelength calibrator. To accurately characterize the parameters of these 
faint lines in the [S II] and [N II] spectra, we followed the same strategy outlined in Haffner (1999) 
for Ha. We first created an average spectrum for each survey "block" (~ 49 pointings taken 
sequentially) . Using this high signal-to-noise spectra, a single Gaussian was fitted to each of the 
atmospheric lines and then subtracted from the individual pointings. Such corrections account for 
hourly variations in these faint lines (typically ~0.02 R). A list of the atmospheric lines in our 
[S II] and [N II] observations is given in Table |l|. The potential atmospheric line at —40 km s~^ in 
[N II] is coincident with interstellar emission from the Perseus arm in all spectra presented here. 
The implications of such a line are discussed in §^. 



3. Results 
3.1. Overview 

Figure || shows samples of Ha, N II, and S II spectra from the region under study. The 
Local (centered near km s~^) and Perseus arm (centered near —40 to —60 km s~^) emission are 
shown as well-separated components in the right-hand panels, where i < 150°, while Local gas 
dominates in the left panels {£ > 150°). Figures §a and |3|a display Ha intensity maps integrated 
over ulsr = — 10 to +10 km s~^ (Figure |2|) and ulsr = —30 to —50 km s~^ (Figure ^), sampling 
the Local and Perseus arms, respectively. The boxes overplotted in them denote the 10° x 12° 
region previously mapped in Ha by Reynolds (1980) and compared to H I observations in detail 
by Reynolds, et al. (1995). The [S II] and [N II] maps look qualitatively similar to Ha and are not 
reproduced here. The line ratios images [N II] /Ha (Figures |2|b and ^j), [S II] /Ha (Figures |2|c and 
^c), and [S II] /[N II] (Figures ^ and ^) are discussed below in §3.3. 
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Two large Ha emitting regions dominate the emission at lower latitudes, above b = —20°. 
The brightest emission arises from NGC 1499 (California Nebula) excited by the 07.51 star ^ Per 
located 540 Iffo pc away {Hipparcos; Ferryman 1997) at i = 160?37, b = — 13?11. Emission in 
the nebula itself, a 5° x 2° crescent centered near £ = 160°, b = —13°, peaks at over 250 R (1 
R = 10^/47r ph cm~^ s~^ sr~^ = 2.4 x 10~^ erg cm^^ s^^ sr^^ at Ha) and is centered around 
■'^LSR = ~5 km s^^, while faint, extended emission above 5 R centered near ^lsr = +5 km 
continues 5-7° away from the nebula. This extended H II region has been studied in some detail 
by Reynolds (1988), who first hypothesized that the bar of emission (/ho ~ 8 R) running from 
e = 150°, b = -20° to £ = 142°, 6 = -10° centered between vlsr = -5 to +5 km s"^ may also be 
ionized by the star. 

The second large emission region, centered near i = 131°, b = —11°, is a superposition 
of separate but equally bright emission regions in the Local and Perseus arms (see Figure [l|) . 
Emission from the Local arm is near vlsr = and has a more circular appearance 

reminiscent of a classical Stromgren sphere (see Figure ^) . This faint H II region has a relatively 
constant brightness of 5-6 R over most of its 13° diameter area. The most promising candidate 
for an ionizing source for this region is (p Per, a B0.5e + sdO binary system located at i = 131?3, 
b = — 11?3 and 220 ^'^32 pc away {Hipparcos; Ferryman 1997; Gies et al. 1998). At that distance, 
the region is 50 pc in diameter. Some details of this region can be found in Haffner (1999). A 
future study (Haffner, Reynolds, & Tufte 1999) examining the properties of B-star H II regions 
in the WHAM Ha Sky Survey will expand on this work. Emission from the Perseus arm in this 
direction at —60 < wlsr < —40 is much more amorphous in shape (Figure ^). 

At higher latitudes and at velocities farther from the LSR, the Ha emission becomes quite 
filamentary. Some of the more notable structures include a faint bar in Figure ^a running from 
i = 150°, b = -33° to ^ = 160°, b = -25°; several, nearly vertical filaments near i = 130°, 
b = —20°; and the shell-like region and associated dissecting filament (more visible at higher 
velocities; see Haffner 1999) centered on the box in Figure The last of these was previously 
studied in detail by Ogden & Reynolds (1985) and Reynolds et al. (1995). 



3.2. Scale Height of the WIM 

Since the Perseus arm gas is well separated in velocity from the Local emission, we can 
explore the decrease of this emission as a function of the height above the Galactic plane. To 
extract information about the scale height from the intensity distribution, we first assume that 
the gas density is primarily a function of z and has the form 

ne{z) = n° e-l^l/-^ cm'^, (1) 
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where is the mid-plane density and H is the scale-hight of the WIM. The Ha intensity is related 
to the gas density by the equation 



2.75Tr/Ha = J ^nidl, (2) 

where (j) is the filling fraction of the emitting gas, T4 is the temperature of the gas in 10^ K, and 
Ina is in R. As long as the temperature, filling factor, and pathlength through the emitting gas, 
J dl = L, are not functions of z we then have 

_ (Pjn'^^fL 2\z\/H _ rO -2\z\/H (o) 

2 75X1^ ~ -'Hq e , [.-i) 

where /j^^ is the intensity at the mid-plane. Finally, substituting z = D tan where D is the 
distance to the arm and b is the Galactic latitude of an observation, and taking the natural log, 
we arrive at 

lnlHa = ln/^„-— tan|6|. (4) 
li 

Figure ^ plots this relationship for a selected region of the map. Each point represents 
the median emission between i = 125° and 152° for each Galactic latitude. The median of the 
emission in each latitude slice is used instead of an average to reduce the effect of local emission 
enhancements {e.g. faint H II regions). The limited longitude range is chosen to avoid the gradual 
blending of the Perseus arm component with the Local gas as the longitude approaches 180°. 
Three radial velocity intervals are presented in the plot, reflecting total (f^LSR = —100 to -1-100 
km s"-*^). Local (vlsr = —25 to -1-100 km s^^), and Perseus arm (ulsr = —100 to —25 km s~-^) 
emission. The integration ranges are expanded from the limited ranges in Figures I and I to 
include the full emission profile in our calculations. The error bars for the Perseus arm are an 
estimate of the physical scatter along a latitude slice calculated from the average deviation about 
the median of emission within that slice. 

Also displayed in Figure ^ is a linear fit to Perseus arm points that are above b = —15°. Closer 
to the plane, H II regions and extinction play a significant role and these points are excluded from 
the fit. The resulting fit parameters are In/y^ = 1.73 it 0.19 and 2D/H = 4.91 =b 0.39, leading 
to a mid-plane intensity for the WIM of 5.7 it 0.2 R and a scale height of 1.0 it 0.1 kpc if the 
Perseus arm is at D = 2.5 kpc (see Reynolds et al. 1995 and references therein). The dotted line 
in Figure Q presents a manual fit to the lower envelope of the Perseus arm emission. The resulting 
parameters are essentially the same: H = 1.0 kpc and = 4.2 R. The 1.0 kpc scale height 
appears to be valid for \z\ from 750 pc to at least 1700 pc, the boundary of this studied region. 

As mentioned above, this estimate assumes that the filling fraction and temperature do 
not depend on If in fact (j) changes as a function of height above the plane, it affects our 
estimate of the electron scale height. As shown below in these new observations suggest that 
the temperature in the WIM rises with distance from the Galactic plane. Such an effect would 
increase our estimate of the scale height. Finally, we also note here that the full role of interstellar 
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extinction has not been taken into account within our restricted fitting region. Future observations 
of H/3 wiU ahow a straightforward correction of the Ha data for this attenuation. A simphstic 
estimate of the magnitude of this correction within our fitted region suggests a decrease in the 
inferred scale height of ~ 10%, somewhat offsetting the expected effect of T(|2;|). 



3.3. Line Ratios 



Figures || and ^ also present pseudo-color images of (b) [N II] /Ha, (c) [S II] /Ha, and (d) 
[S II]/[N II] line ratios integrated over vlsr = —10 to +10 km s~^ (Figure |2|) and i'lsr = —30 to 
—50 km s^^ (Figure sampling the Local and Perseus arms, respectively. The [N Il]/Ha and 
[S II] /Ha maps in both arms generally anti-correlate with Ha emission except for the (p Per H II 



region centered near i = 131°, b = —11° mentioned in §3J above. [S 11] /[N II] ratios are generally 
lower in the the two brighter H II regions in the Local arm, but otherwise show little additional 
correspondence to Ha maps. 

In Figures |5| and [S II] /Ha and [N II] /Ha are plotted against Ha intensity. The radial 
velocity ranges for each panel of the figures are restricted to separate the Local and Perseus arm 
emission (see Figure |l|). As above, the longitudinal extent of the Perseus observations is limited 
to i < 152°. We denote the pointings within the newly discovered (p Per H II region (see ^3.1| ) 
with a square symbol. Note that these points appear to form the upper envelope of the [S II] /Ha 
ratios for their Ha intensities (Figure ||) and actually form a distinctly different horizontal branch 
in the [N II] /Ha ratios (Figure ^). In general, both figures show a marked increase in these ratios 
with decreasing Ha intensity. For comparison to a more traditional H II region in this mapped 
area, the horizontal, solid line at the right edge of both figures shows the average ratio from the 
six pointings with lua > 50 R near NGC 1499: [S Il]/Ha = 0.14 and [N Il]/Ha = 0.39. Ratios for 
a brighter H II region, our calibration source NGC 7000 {lua = 800 R), are [S II] /Ha = 0.06 and 
[N Il]/Ha = 0.20 (see also Reynolds 1985a). 

Figure |^ plots these two line ratios against one another. In both the Local and Perseus arms, 
the ratios appear to be strongly correlated. In this figure, the (j) Per (B-star) and Per (0-star) 
H II regions are highlighted with square and triangle symbols, respectively. Both appear to occupy 
specific loci that are different not only from one another, but from the general WIM background. 
The Perseus arm panel is again restricted to ^ < 152°. 



As discussed in §3.2, the general trend in Ha intensity is that it decreases as the observed 
gas lies farther from the Galactic plane. Figure |8| shows that the relationship of the [S II] /Ha 
and [N II] Ha ratios versus Ha intensity described above translate directly into an increase of 
these ratios with increasing \b\ (or \z\) in the Perseus arm. The figure shows the three line ratios 
[S II] /Ha, [N II] /Ha, and [S II] /[N II] as a function of tan [6|. Both panels are restricted to the 
longitude range i = 125°-152° for consistency. Note also that the [S II]/[N II] ratio is relatively 
constant, with perhaps a slight increase in the Local arm with increasing \b\. 
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4. Discussion 

The tight correlation of [S II] and [N II] emission in Figure ^ results in a fairly constant 
value of [S II] /[N II] over a wide variety of Ha intensities and latitudes in both arms. This result 
has been seen in numerous studies of the extragalactic WIM analogue, commonly referred to as 
the diffuse ionized gas (DIG) (Otte & Dettmar 1999; Rand 1998; Wang, Heckman, & Lehnert 
1997; Greenawalt, Walterbos, & Braun 1997). As Rand covers in detail, current photoionization 
models in which variations in [S II] /Ha and [N II] /Ha arise primarily from variations in the 
ionization parameter {e.g. Domgorgen & Mathis 1994; Sokolowski 1994; see also Bland-Hawthorn, 
Freeman, &: Quinn 1997) have difficulty keeping this ratio constant as [S II] /Ha and [N II] /Ha rise. 
Abundance gradients and H II region contamination can be a larger problem in interpreting such 
ratios in extragalactic studies with the large physical sizes represented by a typical data point. 
Our new observations within the Galaxy provide smaller scales and adequate velocity resolution 
to separate radial components of gas, avoiding some of these problems. 

We suggest that the tight correlation of [S II] /[N II] seen in our data arises simply from the 
dominance of S"*" and ionization states in this gas together with the similar excitation potential 
of these two emission lines. With increasing electron temperature, [S II] /Ha and [N II] /Ha ratios 
rise together with no change in the [S II] /[N II] ratio. Using this premise, namely, that variations 
in these ratios are due primarily to variations in electron temperature, we attempt to derive 
estimates of the electron temperature and ionization state of S within our sample. We start with 
the basic equation for the intensity of emission lines from collisionally excited ions (Osterbrock 
1989), 

. / . 1 2 u fv f 8.63 X 10-*^ S7(i,i) „ 

h (ph s"^ cm-2 ster-i) = ^ j mrie ^ e" w dl, (5) 



where is the collision strength of the transition, uJi is the statistical weight of the ground 

level, Eij is is the energy of the upper level of the transition above the lower, and fi, is the fraction 
of downward transitions that produce the emission line in question. We adopt the following 
parameterized representations of the collision strengths supplied by Aller (1984), valid from 5000 
K to 20,000 K: 

S+:0(%/2,2 1)5/2) = 4.19 X r^-o-osa (6) 
N+ : n{^P,^ D) = 2.68 x T^-^'^^ (7) 

where T4 is in units of lO'' K. A term to account for extinction has not been included in Equation]^, 
since the three lines we are comparing are close in wavelength. 

To arrive at any estimates, we must be able to compute the path averaged integrals in these 
equations. Since there is little knowledge yet of the details of the spatial distribution of the 
physical parameters in the WIM, we take the simplest approach and assume that physical values 
do not vary along a path. For the ratios we calculate below, we need to make this assumption for 
the temperature, ionization, and abundance of the gas. Our estimates should then be viewed as 
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representative values, rather than accurately reflecting local conditions from where the emission 
arises. Since we are able to cleanly separate Local and Perseus arm gas, large-scale Galactic 
variations can be probed along lines of sight, even with this assumption. This simplification 
also implies that the line emission from the three ions originates from the same location. This 
is probably a valid assumption in the WIM, unless local processes on small scales {e.g. shocks) 
dominate the ionization structure. In the brighter, discrete H II regions, however, this assumption 
most likely breaks down since the ionization structure of the nebula could be changing significantly 
within our beam. 

With this simplification, we arrive at these two equations for the intensity of [S II] and [N II] 
emission, 

76716 (R) = 2.79 X 10^ ' [I] ] T^-O-^^^ EM, (8) 



and 



H+\ ^ ( N 



/6583 (R) = 5.95 X loM ^ (9) 



where T4 is in units of 10^ K and EM = J nidi is the emission measure in units of cm ^ pc. Note 
also that / is now in R. The [S II] /[N II] intensity ratio is then 

-^6716 ^ gg {h) (~) g0.04/T4 2^-0.119 _ ^-^Q^| 



'6583 



f N\ ( N+ 
\h) \ N 



We use the solar (S/H) = 1.86 x 10"^ from Anders & Grevesse (1989) and (N/H) = 7.5 x 10"^ 
from Meyer, Cardelli, & Sofia (1997) for the gas-phase abundances of S and N. 

Because of the similar first ionization potentials of N and H (14.5 and 13.6 eV) and a weak 
charge-exchange reaction, N+/N'^ tracks H+/H*^ in photoionization models (Sokolowski 1994), so 
that in Ha emitting regions, where the fractional ionization of H has been measured to be near 
unity (Reynolds et al. 1998), — > 1 also. However, with a second ionization potential of 29.6 
eV, N most likely ionizes no higher than N^. Recent photoionization models run by Howk & 
Savage (1999) over a wide range of temperatures and ionization parameters show that even if the 
ionizing spectrum in the WIM is as hard as Tg ~ 40, 000 K, < 0.3. Furthermore, Reynolds 

& Tufte (1995) and Tufte (1997) found that He+/He ranges from 0.3 to about 0.6 in the WIM, 
suggesting a spectrum that is significantly softer than 40,000 K, and therefore implying a very 
small N"'"^/N ratio. Sulfur, on the other hand, has a low enough first ionization potential, 10.4 
eV, that it is most likely fully ionized in the gas we are studying. In addition, its second potential, 
23.4 eV, lies slightly below the neutral He edge at 24.6 eV. The existence of He"'" in the WIM 
(Reynolds & Tufte 1995; Tufte 1997) and the model calculations of Howk &; Savage (1999) both 
imply that indeed some S will be ionized to . 

In the limit where S"'"/S and N^/N are both unity and at a typical WIM temperature of 



T4 = 0.8, Equation 10 gives /67i6/-^6583 = 1-26. This relationship is plotted as the solid line in 
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Figure ^. Dotted lines of lower slope show the effect of decreasing S+/S (presumably due to an 
increase in S^^) while keeping N+/N = 1. In this interpretation, the data imply that S"^/S ranges 
between about 0.25 and 0.8, with an average value for the Local diffuse background of 0.6-0.65. 
Note that the location of the H II region symbols (filled triangles and squares) in Figure ^ indicate 
that these regions have systematically lower S^/S than the diffuse background. In these cases the 
lower ratios are more likely due to the expected behavior that closer to the ionizing source, where 
the ionization parameter is large, more S"*" is ionized to producing a corresponding decrease 
in the observed [S II] intensity relative to Ha (Domgorgen &; Mathis 1994; Sokolowski 1994). 

Just as the [N II] /Ha and [S II] /Ha ratios can together give us information on the ionization 
state of S, examining the ratios alone can provide some estimate of the temperature. From 
Equations |^ and ^ the [N II] /Ha ratio is given by 




= 1.63 X loM ^ 1 I ^ 1 I ^ 1 e~'-''/^^. (11) 



Since N+/N ~ H+/H (see above), the H and N ionization factors can be removed from Equation 11 
with the result that for a a given (N/H) abundance, [N II] /Ha is only a function of the electron 
temperature. As an example, the typical WIM temperature of r4 = 0.8 then leads to [N II] /Ha 
= 0.72. Lines of constant temperature are plotted on Figure ^ as dashed vertical lines. With a 
better estimate for temperature from this line ratio, we also recompute Equation |l^, including the 
temperature dependencies this time, to form the dash-dotted line in Figure |9[ This correction is 
small. 

Figure ^ shows the same grid of S''"/S and T overplotted on the Perseus arm ratios. The 
limits placed on S''"*"/S in the WIM by Howk & Savage (1999) are at odds with the significant 
number of points below the S^/S = 0.5 line in the Perseus arm. As noted in there may be an 
atmospheric line within the integration range of Perseus arm [N II] emission. By comparing [S II] 
and [N II] spectra toward the same locations, we derive an upper limit of 0.1 R for the intensity 
of this line. The maximum effect of removing this line is displayed for a few sample points in 



Figure 10 by translating the solid sample error bars points to the location of the dotted error 
bars. Such a correction would imply a slightly higher S^/S ratio on average for the Perseus arm, 
but not enough to explain the entire difference between the arms. We may be probing a harder 
radiation field in the Perseus arm gas than their model input spectra or our assumed value of S/H 
may be too high. A reduction in the gas phase sulfur abundance by 0.5-0.75 (0.3-0.1 dex) would 
bring our inferred S'''^/S fractions more in line with their models. 

The Perseus arm gas appears to have systematically lower S^/S than the Local arm, with 



an average value closer to 0.5. Equation |10| shows that differences in abundances between the 
Local and Perseus arm could change the slope of the [S II] /Ha and [N II] /Ha relationship. 
However, Afflerbach, Churchwell, & Werner (1997) measured the Galactic abundance gradients in 
ultra-compact H II regions for O, S, and N. Using their best fits, N/H falls off faster than S/H 
with increasing Galactocentric distance so that S/N is ~ 5% higher in the Perseus arm. This is 
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the opposite effect needed to explain tlie difference in tlie Perseus and Local arm emission ratios. 
More likely, Figure ^ suggests that this region of the Perseus arm has a higher fraction of S"*""*" 
than the Local arm region. Rand (1998) observed an increase in [O IIl]/H/3 with height above 
the plane of NGC 891 and concluded that an additional source of ionization is needed beyond 
photoionization from hot stars in the disk. In a region where O^^ ions are being maintained, S"*" 
becomes a less dominant species. This is consistent with the inferred differences in S^/S between 



the Perseus arm (Figure 10) and Local arm (Figure ^ data, since lines of sight pass through 
the Perseus arm at much higher \z\ than for Local gas. One problem with this interpretation is 
that we do not see a substantial fall in the [S n]/[N II] ratio with increasing \z\ in the Perseus 
arm (Figure |8|). Instead of an additional ionizing source, it is also possible that there are marked 
differences in the strength of the halo ionizing field between the Local and Perseus arms. For 
example, NGC 891 shows considerable variation in the extent of extended Ha emission along the 
galactic plane (Rand, Kulkarni, & Hester 1990; Rand 1998). 

The scatter of data points suggests significant variations in the electron temperature within 
the mapped region, from about 6,000 K to 10,000 K. The most striking result from these 
temperature estimates comes when re-examining the trend of the [N II] /Ha ratio in Figures |^ 
and 1^. From the rise in this ratio with increasing \z\ and Equation |ll| (or the grid of Figure |l0|), 
we infer that the temperature in the halo rises from about 7,000 K at \z\ = 0.75 kpc {b ~ —17°) 
to over 10,000 K at \z\ = 1.75 kpc (6 ~ —35°). If the electron temperature is truly rising with \z\ 



as our results suggest, the calculation of the scale height of Ha emitting gas (§|3.2|) becomes more 



complicated. If we confirm this effect with future observations (see the net effect (Equation |2|) 



is that the true scale height of the gas is larger than the value calculated in §3.2 by about 20%. 



Note also that since fainter regions tend to have higher [N II] /Ha (Figure ^), our analysis 
implies that fainter regions are hotter. Since EM oc /ho^'^'^) higher temperature regions naturally 
produce lower Ha intensities for a given EM; however, even after converting /h« to EM using 
a varying temperature estimate from the [N II] /Ha observations, regions with lower EM are 
still hotter. Since EM depends on ^rig and the pathlength through the emitting region, one of 
the two must be smaller in these fainter, hotter regions. Pathlength is an inherently unphysical 
parameter to link to temperature and we can think of no reason why our data set would bias 
shorter pathlength sightlines toward higher temperature regions. Thus, we conclude that the 
hotter regions most likely have smaller {n^)- Since Ue in our mapped region of the Perseus arm 
decreases with increasing \z\ (see ^3.2| ), we are not able to ascertain from these data whether rig or 
\z\\s the fundamental parameter influencing the electron temperature. 

One possible explanation for the discrepancy between these results and photoionization 
models {e.g. Domgorgen &: Mathis 1994; Sokolowski 1994) is that at high \z\ or low Ue an 
additional heating source may dominate over the heating from photoionization (Bland-Hawthorn 
et al. 1997). Since the heating per volume from photoionization is proportional to rig (it is limited 
by recombination), a heating term that is proportional to or did not depend on rig at all would 
become dominant at sufficiently low densities. Such an effect would decouple the heating from the 
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ionization of atoms, keeping [S 11] /[N II] constant, while driving up the [S II] and [N II] intensities 
relative to Ha. We also note that if such a process exists in the halo, the increase in [O III]/H/3 
seen by Rand (1998) in the halo of NGC 891 could be a direct consequence of this additional 
heating source, and would not require a secondary ionization source. Possible heating sources 
include, for example, the dissipation of interstellar turbulence (Minter & Spangler 1992) and the 
photoionization of interstellar grains (Reynolds & Cox 1992). The idea of an additional heating 
source will be explored further in a related work (Reynolds, Haffner, & Tufte 1999). 

5. Summary 

We have presented the first comprehensive, velocity-resolved maps of [N II] and [S II] in the 
Galaxy. Combining the Ha maps from the WHAM Ha Sky Survey in this region, we draw the 
following conclusions: 

1 . We confirm earlier estimates of the distribution of Ha emission in the Perseus arm (Reynolds 
1997). With continuous coverage from b = —6° to —35° {\z\ = 300 to 1750 pc at d = 2.5 
pc), we find that the electron distribution derived from the emission is described well by an 
exponential with a scale height of H = 1.0 ± 0.1 kpc. 

2. [S II] /Ha and [N II] /Ha emission in our Galaxy shows the same rise with increasing \z\ 
as found in observations of edge-on galaxies (Rand 1998; Golla et al. 1996). These ratios 
appear to be correlated with each other over a wide range of emission measures and \z\. 

3. We suggest that the correlation between [N II] /Ha and [S II] /Ha can be attributed to the fact 
that changes in these ratios are primarily due to changes in electron temperature rather than 
ionization parameter, and that scatter in the relationship is caused by variations in S^/S. 
The maps of [N II] /Ha and [S II] /[N II] (Figures § and ^) can then be viewed as distributions 
of electron temperature and ionization state, respectively. With this interpretation, we 
derive temperatures in the WIM of 6,000 K to 10,000 K and S+/S ratios of 0.3 to 0.8. The 
average value of S^/S appears to be lower in the Perseus arm. Higher ionization and lower 
temperatures are mostly found in the two H II regions surrounding the O star ^ Per and the 
B star system cp Per. Differences in temperature between regions may be able to be directly 
confirmed with observations of the faint, higher-level [N II] A5755 transition, which when 
combined with our A6483 measurements can directly measure the electron temperatures in 
the emitting gas. 

We thank Kurt Jaehnig and Jeff Percival of the University of Wisconsin's Space Astronomy 
Lab for their dedicated engineering support of WHAM; Nikki Hansen, Mark Quigley, and Brian 
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Fig. 1. — Sample Ha (solid line), [N II] (dashed line), and [S II] (dotted line) spectra from several 
directions are plotted as intensity versus LSR velocity. The total integrated intensity of the emission 
is labeled in the upper-left corner of each plot. The dashed vertical line denotes the delineation for 
integration ranges between the Local and Perseus arms in this paper. 

Fig. 2. — (a) Ha emission from the local Orion Arm. This pseudo-color image shows the Ha 
emission integrated over ulsr = —10 to +10 km s~^. The axes are Galactic longitude and 
latitude. The overplotted box shows the location of the original Ha background survey in this 
region (see text). Emission line ratio maps of (6) [N Il]/Ha, (c) [S Il]/Ha, and (d) [S n]/[N II] 
are also shown for the same velocity band. The color scale in each image is histogram equalized to 
enhance spatial detail. {Preprint Note: A color PostScript version of this figure can be found at 
http : //www . astro . wise . edu/wham/ papers . html) 

Fig. 3. — (a) Ha emission from the Perseus Arm, approximately 2.5 kpc distant. This pseudo-color 
image shows the Ha emission integrated over fLSR = —50 to —30 km s^^. The axes are Galactic 
longitude and latitude. The overplotted box shows the location of the original Ha background 
survey in this region (see text). Emission line ratio maps of (6) [N Il]/Ha, (c) [S II]/Ha, and (d) 
[S II] /[N II] are also shown for the same velocity band. The color scale in each image is histogram 
equalized to enhance spatial detail. {Preprint Note: A color PostScript version of this figure can 
be found at http://www.astro.wisc.edu/wham/papers.html) 

Fig. 4. — The natural logarithm of the median Ha intensity between i = 125° and 152° is plotted 
versus the tangent of Galactic latitude. Total (vlsr = —100 to +100 km s~^). Local (fLSR = —25 
to +100 km s^^), and Perseus arm (wlsr = —100 to —25 km s^^) emission are denoted by squares, 
diamonds, and vertical bars, respectively. A linear best fit to the Perseus arm emission above 
b = —15° is displayed as a solid line. The dotted line is a fit to the lower envelope of the data (see 

§iD. 

Fig. 5. — [S II] /Ha is plotted versus Ha intensity for both the Local and Perseus arm emission. 
For the Perseus arm, the longitude range has be restricted to ^ = 125° through 152°. Crosses 
show representative error bars for select points. The uncertainty in the ratio at 0.5 R is 7-8%. 
Filled squares denote pointings within 6.5° of the center of the Per H II region mentioned in § |3.1| : 
£ = 131°, b = —11°. The solid line denotes the average ratio of pointings with Ina > 50 R (i.e. 
"classical" H II regions). 

Fig. 6. — Same as Figure H, but for [N II] /Ha versus Ha. 

Fig. 7. — Same as Figure but for [S II] /Ha versus [N II] /Ha. Additional plot symbols (triangles) 
denote pointings within 5° of the (rough) center of the ^ Per H II region: £ = 160°, b = —15°. 

Fig. 8. — Three emission line ratios versus the tangent of Galactic latitude: [S II] /Ha (dotted 
line), [N II] /Ha (dashed line), and [S II] /[N II] (solid line). Each data point represents the ratio of 
the median emission between i = 125° and 152°. Intensities have been integrated over the velocity 
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range ^^lsr = —25 to +100 km s for the top panel and ulsr = —100 to —25 km s in the bottom 
panel. 

Fig. 9. — Same as Figure ^ for the Local arm only, except with a grid for inferred temperatm'es 
and S^/S ratios. 

Fig. 10. — Same as Figure |7| for the Perseus arm only, except with a grid for inferred temperatures 
and S^/S ratios. As before, solid crosses denote error bars for representative samples; dotted 
crosses show how the positioning of these data points would change if a weak terrestrial line at the 
velocity of the Perseus arm is as bright as 0.1 R (see text). 
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Table 1. Atmospheric Lines in [S II] and [N II] Observations 
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+26.8 





0.06-0.38 


+25 





0.05-0.10 
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10 
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-76 


10 


0.00-0.05 


-113.2 
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(-40 





0.05-0.10)^ 



^With respect to the geocentric center of the emission hne. 

^FuU variation in intensity seen throughout the observation period. Hourly variations are ~ 0.02 

R. 

'^Contaminated by Perseus arm emission. See text. 




I I — I — I — I — I — I — I — I — I-Ih — I — I — I — I — I — I — I — I — I 




-100 -50 50 -100 -50 

Velocity [km s"'] Velocity [km s 



- 20 - 



Figure 2 is available separately as a GIF from LANL Preprints or as a color PostScript file from 

htt p : / / www, astro . wise . edu/ wham / papers . html 
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Figure 3 is available separately as a GIF from LANL Preprints or as a color PostScript file from 

htt p : / / www, astro . wise . edu/ wham / papers . html 
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